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ABSTRACT 



Aims. In a previous study we suggested that the broad-band emission and variability properties of BL Lacertae can be accounted for by a double 
synchrotron emission component with related inverse-Compton emission from the jet, plus thermal radiation from the accretion disc. Here we 
investigate the matter with further data extending over a wider energy range. 

Methods. The GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) monitored BL Lacertae in 2008-2009 
at radio, near-IR, and optical frequencies to follow its flux behaviour. During this period, high-energy observations were performed by XMM- 
] Newton, Swift, and Fermi. We analyse these data with particular attention to the calibration of Swift UV data, and apply a helical jet model to 
O ii interpret the source broad-band variability. 

I ■ Results. The GASP-WEBT observations show an optical flare in 2008 February-March, and oscillations of several tenths of mag on a few-day 
O ' time scale afterwards. The radio flux is only mildly variable. The UV data from both XMM-Newton and Swift seem to confirm a UV excess 
I that is likely caused by thermal emission from the accretion disc. The X-ray data from XMM-Newton indicate a strongly concave spectrum, 
CZ2 , as well as moderate (~ 4-7%) flux variability on an hour time scale. The Swift X-ray data reveal fast (interday) flux changes, not correlated 
with those observed at lower energies. We compare the spectral energy distribution (SED) corresponding to the 2008 low-brightness state, which 
was characterised by a synchrotron dominance, to the 1997 outburst state, where the inverse-Compton emission was prevailing. A fit with an 
I \ inhomogeneous helical jet model suggests that two synchrotron components are at work with their self inverse-Compton emission. Most likely, 
they represent the radiation from two distinct emitting regions in the jet. We show that the difference between the source SEDs in 2008 and 1997 
can be explained in terms of pure geometrical variations. The outburst state occurred when the jet-emitting regions were better aligned with the 
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f ^ ' line of sight, producing an increase of the Doppler beaming factor. 
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Conclusions. Our analysis demonstrates that the jet geometry can play an extremely important role in the BL Lacertae flux and spectral variability. 
Indeed, the emitting jet is probably a bent and dynamic structure, and hence changes in the emitting regions viewing angles are likely to happen, 
with strong consequences on the source multiwavelength behaviour. 

■ Key words, galaxies: active - galaxies: BL Lacertae objects: general - galaxies: BL Lacertae objects: individual: BL Lacertae - galaxies: jets 
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1. Introduction iRavasio et ail d2002l); iBottcher et ail d2003l): IVillata et alj 



(2004ha); Bach et all d2006l); iPapadakis et all d2007l) 



V , Blazars are active galactic nuclei whose extreme properties are viUa ^fJ J 2 Q09bl) : lLarionovtttn (2010). 

r »_J thought to be owing to their relativistic lets pointing toward us. T . U .. • . , I \U n nnk 1 A 1 

Q & o , „f, , In a recent paper, Raiteri et al. (2009) analysed the mul- 

. 03 BL Lacertae, the prototype of the BL Lac objects blazar class tiwaV elength data from the 2007-2008 WEBT campaign, in- 

has been the target of many campaigns by the Whole Earth dudi three inti b XMM-Newton. The XMM-Newton 

Blazar Telescope (WEBT) collaboration] since 1999. The tens , . , , mr u- u »j»uj 

K. . , ' ,. , , , data revealed a UV excess, which was interpreted to be due 

of thousands of optical-to-radio data collected by the WEBT tQ emisfjion from ^ accretion ^ afj weU afj a 

allowed us to study its multiwavelength flux variability colour ^ curvature in the x band xhe authors constmcted spec- 

behaviour, and the correlations among flux variations in different . , ■ • . ,. , ccr , s tDI , , 

, , ., , • ,. • r , ,. tral energy distributions (SEDs) of BL Lacertae corresponding 

bands, and revealed a possible periodicit y of the rad i^out- tQ yarious hg where ^ SQurce wag ^ different bri htness 

bursts. The results have b een published by| Villata et al. | d2002j); using both their Qwn data and data from ^ literature 

Send offprint requests to: C. M. Raiteri They applied the inhomogeneous, rotating helical jet model by 

* The radio-to-optical data presented in this paper are stored iMafa & Raiteri ! 1999, see also |Raiteri et al. || l999j | Raiteri et al. | 

in the GASP-WEBT archive; for questions regarding their avail- |2003|, |Ostorero et al.||2004|) to fit the SEDs, and suggested that 

ability, please contact the WEBT President Massimo Villata the broad-band spectral properties of BL Lacertae may result 

(villata@oato . inaf . it). from the combination of two synchrotron emission components 

1 [http : //www . oato . inaf . i t/blazars/webt/l with their self inverse-Compton emission, plus a thermal com- 
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Fig. 1. R-band light curve by 
the GASP collaboration from 
February 2008 to February 
2009. Data are not corrected 
for the host galaxy contamina- 
tion or Galactic extinction. The 
participating observatories are 
marked with different symbols 
and colours. The total number 
of data points is indicated in the 
upper right. 



ponent from the disc. Subsequently, Capetti et al. (2010) anal- 
ysed optical spectra acquired in the same period with the 3.56 
m Telescopio Nazionale Galileo (TNG). They found a broad Ha 
emission line, with luminosity of ~ 4 x 10 41 erg s -1 and FWHM 
of ~ 4600 km s" 1 , even br ight er than that found in 1995 -1997 
bv lVermeulen et al.1 (1 19951) and lCorbett etafl ([1996, 2000). This 
favours the hypothesis that the UV excess is caused by ther- 
mal emission from the accretion disc, the most likely source 
of ionising photons for the broad line region. The multi wave- 
length data available for the|Raiteri et al. (2009) analysis lacked 
simultaneous information in the y-ray band, so that the inverse- 
Compton spectral region was poorly constrained. But in 2008 
the F ermi satellite was able to detect BL Lacertae dAbdo et al.L 
I2010ah . even if in a low state compared to the pas t detections by 
the Compton Gamma Ra y ObservatorjQ (CGRO, Hartma net al.l 
l!999llBToom et al.ll997l) . In the same period, observations in the 
UV and X-ray bands were performed by Swift, while in the op- 
tical, near-IR, mm and cm radio bands the source was monitored 
by the GLAST- AGILE Support Program (GASP) of the WEBT. 
This offered the unique opportunity to study the source emission 
over a very extended spectral range. The results of this further 
investigation effort on BL Lacertae are presented in this paper. 

2. GASP observations 

The GASP was born in 2007 as a WEBT project, with the aim 
of monitoring a list of 28 y-ray loud blazars in the optical, near- 
IR, mm, and cm radio bands during the y-ray observations of 
the AGILBl and Ferrnfl (f ormerly GLAST) satellites (see e.g. 
IVillata et all [2008, 20093)- Data are collected periodically by 
the WEBT President, who checks the consistency of the vari- 
ous datasets. The GASP light curves are then available for mul- 
tiwavelength studies, mostly in the framework of the GASP 
collaboration with the AGILE and Fermi research teams. The 
GASP data presented in this paper were taken at the observato- 
ries listed in TableQ] 

The optical data were calibrated with respect to a com- 
mon choice of refer ence stars in the same field of the sour ce 
dBertaud et al.ll 19691 in U and B bands: iFiorucci & Tostil[l996l in 

2 http : //heasarc . gsf c . nasa . gov/docs/ cgro/ 

3 |http : / / agile . ias f-roma . inaf . It/1 

4 http://fermi.gsfc.nasa.gov/ 



Table 1. List of optical, near-IR, and mm-cm radio observatories 
contributing data to this work. 



Optical and near-infrared 



Observatory Tel. size Bands 

[cm] 

Abastumani, Georgia 70 R 

Armenzano, Italy 35 BR1 

Armenzano, Italy 40 BVRI 

Calar Alto, Spain" 220 R 

Campo Imperatore, Italy 1 10 JHK 

Crimean, Ukraine 70 BVRI 

El Vendrell, Spain 20 R 

Kitt Peak (MDM), USA 130 UBVRI 

L'Ampolla, Spain 36 R 

Lulin, Taiwan 40 R 

New Mexico Skies, USA 30 VRI 

Roque (KVA), Spain 35 R 

Sabadell, Spain 50 R 

St. Petersburg, Russia 40 BVRI 

Talmassons, Italy 35 BVR 

Teide (BRT), Spain 35 BVR 

Tuorla, Finland 103 R 

Valle d'Aosta, Italy 81 BVRI 

Radio 

Observatory Tel. size Frequencies 

[m] [GHz] 

Mauna Kea (SMA), USA Wx(F 230, 345 

Medicina, Italy 32 5, 8, 22 

Metsahovi, Finland 14 37 

Noto, Italy 32 43 

UMRAO, USA 26 4.8, 8.0, 14.5 



" Calar Alto data were acquired as part of the MAPCAT (Monitoring 
AGN with Polarimetry at the Calar Alto Telescopes) project. 
* Radio interferometer including 8 dishes of 6 m size. 



V, R, and /). The source photometry was evaluated from a circu- 
lar region with an 8 arcsec aperture radius, while the background 
was taken in a surrounding annulus with 10 and 16 arcsec radii. 
In this way the measure is essentially seeing-independent and all 
datasets are affecte d by the same cont amination from the light 
of the host galaxy. Rait eri et~afl (|2009) estimated that with the 
above prescriptions the contamination amounts to about 60% 
of the host total flux density, which is 0.36, 1.30, 2.89, 4.23, 
5.90, 11.83, 13.97, and 10.62 mJy in the U, B, V, R, I, J, H, 
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Fig. 2. June-November 2008 light curves of BL 
Lacertae from UV to near-IR after correction 
for the host galaxy contribution, but not for the 
Galactic extinction. Data taken by the GASP- 
WEBT collaboration are plotted as blue circles, 
UVOT data as red triangles. The UVOT u, b, 
and v light curves have been shifted to match 
the ground-based U, B, and V ones (see text for 
details). 



and K bands, respectively. When converting magnitudes into 
flux densi ties, we corrected f or the Galactic extinction accord- 
ing to the lCardellietail (1 19891) laws, using R y = 3.1, the stan- 
dard value for the d i ffuse i nterstellar medium, and Ag = 1 .42 
(from|Schlegel et all Il998l) . We adopted the absolute fluxes by 
iBessell et al.ldl998h . 

Figure Q] shows the best-sampled total /?-band light curve 
from February 2008 to February 2009 built with GASP data, 
which are not corrected for the host galaxy contribution here. 
A noticeable flare was observed at the beginning of the period, 
in 2008 February-March; afterwards both the average bright- 
ness level and the variability amplitude decreased. However, 
the source remained active, its brightness oscillating by several 
tenths of magnitude on a few-day time scal e. This is not an un- 
usual behaviour for BL Lacertae (see e.g. Raiteri et al., 2009, 
who reported on a 0.9 mag brightening in 24 hours). 

Optical data at other wavelengths as well as near-IR data are 
shown in Fig. [2] for the period June-November 2008 (the UV 



data displayed in the figure are presented in Sect. 0. The host 
galaxy contribution has been subtracted to distinguish the be- 
haviour of the active nucleus. This reveals that the brightness 
evolution follows the same trend in the various bands, but mag- 
nitude variations are more pronounced at higher frequencies, 
which is a common feature of BL Lac objects. For example, the 
brightness increase following the almost symmetric dip around 
JD ~ 2454681 involved a variation of 0.79, 0.64, 0.57, and 0.54 
mag in ~ 5 days in the B, V, R, and / bands, respectively. We 
notice that the source redshift is z = 0.0686 dVermeulen et al.L 
1995) and hence the broad Ha emission line enters the tails 
of the R and / passbands. However, referring to Capett iet al.l 
d2oToh . one can estimate that its flux contribution is only a few 
thousandths of that of the continuum. Hence, the presence of the 
line cannot affect the variability in these bands, at least when the 
source brightness is at these levels. Another interesting exam- 
ple of fast variability is the rise of ~ 0.9 mag in 4 days, from 
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JD = 2454752.3 to 2454756.3, in the R band, which unfortu- 
nately was not observed in other bands. 

Radio flux densities at cm-mm wavelengths are displayed 
in Fig. [3] together with de-reddened and host-galaxy subtracted 
/?-band optical flux densities. The former are complemented by 
data from the VLA/VLBA Polarization Calibration Database^ 
(PCD). The flux variation amplitude appears to decrease from 
the highest to the lowest frequencies, as usual. One interest- 
ing feature is the fast radio flare that is visible in the 37 GHz 
light curve at JD ~ 2454760, because it occurred simultane- 
ously with an optical flare. The correlation between optical and 
radio flux variations in BL Lacertae has been the subject of sev- 
eral studies, and most of them found a correlation with a long 
time delay (a fe w months) of the radio after the optical flux 
changes ("see e.g. | Hufnagel & Bregmarlll992l;lTornikoski et all 
1994 b] IClements et all Il995b IVillata et all l2004at iBach et all 
2006; [Villata et all 2009b). However, sim ultaneous variations 



have already been found ( Tornikoski et all[l994al) . The 37 GHz 
data we are dealing with are affected by large uncertainties due 
to unfavourable weather conditions; but an increase of the radio 
flux is visible also at 43, 22, and 14.5 GHz, giving strength to 
the possibility that these events are correlated. We also notice 
that there is neither a contemporaneous nor a delayed clear radio 
counterpart to the optical flare observed at the beginning of the 
period (JD ~ 2454520). Yet, one would have expected to see it 
in the high-frequency radio light curves. A better sampling in the 
mm wavelength range perhaps would have helped to understand 
whether the variability mechanism responsible for this flare af- 
fected the optical emission region only or if it extended also to 
the radio domain. 

3. Swift observations 

The Swift satellite observed BL Lacertae in 2008 August, 
September, and October, for a total of 24 pointings. In particular, 
from August 20 to September 9 a daily sampling was obtained. 

3.1. UVOTdata 



The Ultraviolet/Optical Telescope (UVOT; iRoming et al.ll2005h 
onboard the Swift spacecraft acquires data in the optical v, b, 
and u bands , as we ll as in the UV filters uvwl, uvm2, and uvw2 
dPoole et alll2008l) . 

We reduced the BL Lacertae data with the HEAsoft pack- 
age version 6.7, with CALDB updated at the end of November 
2009. Source counts were extracted from a circular region with 
a radius of 5 arcsec, while background counts were estimated in 
a neighbouring source-free region. When multiple exposures in 
the same filter were present during an observing epoch, we first 
processed each frame separately with the task uvotmaghist and 
then binned the results. These values were compared with those 
obtained by first summing the frames acquired in the same band 
with uvotimsum, and then performing the aperture photometry 
with the task uvotsource. We verified that the two methods are 
equivalent. 

The final UVOT light curves are shown in Fig. [2] We sub- 
tracted the host galaxy contribution, taking into account that 
with the 5 arcsec aperture radius we used for the photometry, 
about 50% of the host flux w as included. We ad opted the host 
galaxy flux densities given bv lRaiteri et alj d2009l see also Sect. 
[2]); these authors also discussed that the host contribution can be 
considered negligible in the UV. 
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Fig. 3. Optical flux densities (R band, top panel), after correc- 
tion for Galactic extinction and host galaxy contamination, com- 
pared to radio flux densities at different frequencies. Blue cir- 
cles represent GASP data; red diamonds indicate data from the 
VLA/VLBA PCD. The vertical green line corresponds to the 
XMM-Newton pointing of May 16-17; the yellow strip high- 
lights the period of Swift observations, from 2008 August 20 to 
October 2. 



The comparison between the UVOT u, b, and v data and the 
U, B, and V data taken by the GASP observers reveals that an 
offset is present between the space light curves and the ground- 
based ones. We estimated mean offsets U — u — 0.2, B - b — 0.1, 
and V — v = -0.05. The UVOT light curves shown in Fig. [2] 
have been shifted accordingly. Taking into account that the av- 
erage UVOT colour indices of BL Lacertae are u - b ~ -0.4 
an d b — v ~ . 8, the above offsets disagree with those derived 
bv lPoole et alj d2008l) for the objects on which they based their 
photometric calibration of UVOT, i.e. Pickles stars and GRB 
models. I ndeed, these object s have a different spectral shape, 
so that the lPoole et alj J2008) calibrations may not hold for BL 
Lacertae. 

The UVOT data confirm the variability trend traced by the 
ground-based ones, extending it to UV frequencies. This in- 
dicates that the variability mechanism affecting the near-IR- 
optical emission, which is dominated by beamed synchrotron 
radiation, can also produce flux changes in the UV, where a con- 
tribution from the synchrotron emission is thus expected, besides 
a possible contribution from thermal disc radiation. 



4 



C. M. Raiteri et al.: Another look at the BL Lacertae flux and spectral variability 



-14.2 
— 14.4 


w2 


m.2 


w1 


u 


b 


v \ 


-14.6 
-14.8 












- 


-15.0 














-15.2 
-15.4 












• <OM> — 

♦ <UV0T> 



3.5 3.6 
log A [A] 



B -10.4 ; 



fc. -10.5 r 



It 



• <0M> 

+ <UV0T> 



14.7 14.8 14.9 15.0 15.1 15.2 

log v [Hz] 

Fig. 4. Top: Observed spectrum of BL Lacertae in the optical- 
UV band. Red diamonds refer to the mean UVOT spectrum re- 
sulting from 16 epochs analysed in this paper. Blue circles rep- 
resent the average spectrum obtained from the three observa- 
tions of the OM instrum ent onboard XMM-Newton in 2007- 
2008 dRaiteri et all 120091) . normalised to the mean UVOT spec- 
trum in the v band. The solid line is the log-parabolic fit used 
in the calibration procedure. Bottom: Optical-UV SEDs of BL 
Lacertae. Blue circles and red diamonds are derived from the 
OM and UVOT average spectra shown in the top panel by using 
standard prescriptions to obtain dereddened flux densities. Black 
squares represent the mean UVOT SED after recalibration as ex- 
plained in the text. In bo th panels the filter labels are centred on 
the standard UVOT A eff dPoole et alll2008h . 



The derivation of the source intrinsic flux densities for fur- 
ther analysis (see Sect. [5]) requires some att ention. In their pa- 
per on the photometric calibration of UVOT [Poole et al] (12008) 
give effective wavelengths of 5402, 4329, 3501, 2634, 2231, and 
2030 A for the v, b, u, uvwl, uvml, and uvwl filters, respec- 
tively, but warn that the A e s of the UV filters will be longer for 
very red spectra. Moreover, they provide count-rate-to-flux con- 
version factors for both Pickles stars and GRB models, but in 
the UV bands their validity range is limited to b - v = 0.1 and 
b - v — 0.03, respectively, while BL Lacertae has b - v ~ 0.8. 

In order to attenuate possible calibration problems, we thus 
calculated both the effective wavelengths A^ and count-rate- 
to-flux conversion factors CF for the UVOT filters by fold- 
ing the BL Lacer tae spectrum with their effective areas (see 
Pool e et all |2008). We first built a composite observed spec- 
trum of BL Lacertae by combining a mean OM spectrum (ob- 
tained from the th ree XMM-Newton pointings of 2007-2008, 
Raite ri et"a"D 12009) with an average UVOT spectrum (resulting 



from 16 UVOT observing epochs analysed in this paper). To 
compensate for the different brightness state, we increased the 
OM flux densities by ~ 6% so that the two spectra match in 
the V band. The composite spectrum is shown in Fig. [4] (top 
panel), together with its log-parabolic fit that we used in the 
foldi ng procedure. Th e resulting effective wavelengths (see Eq. 
8 in lPoole et al.1 120081) are: 5439, 4381, 3500, 2776, 2295, and 
2225 A for the v, b, u, uvwl, uvml, and uvwl filters, respectively, 
showing a clear shift towards longer wavelengths in the ultravio- 
let. As for the count-rate-to-flux conversion factors, we obtained 
2.60, 1.47, 1.65, 4.31, 8.54, and 6.72 x 10~ 16 erg cirT 2 s" 1 A" 1 
from the v to the uvwl band, respectively . These new CF dif- 
fer from those given bv lPoole et al.1 d2008) for the GRB models 
by ^ 1%, with the only exceptions of CF(uvwl) and CF(uvw2), 
which are now 8% largefl 

The new A e g would produce a decrease of extin ction in the 
uvwl and uvml bands with respect to those given bv lPoole et al.l 
d2008h . and an increase in the uvwl band. Indeed, the Galactic 
mean extinction curve shows a dramatic bump peaking at A ~ 
2175 A owing to absorption by graphite dust. Actually, an ac- 
curate evaluatio n of extinction in thi s critical frequency range 
requires that the lCardelli et all (1 19891) law is folded with the fil- 
ter's effective area and BL Lacertae spectrum, similarly to what 
was done above for the A e $ and CF: 



A a =2.5 log 



J dAE^FAA)!^' 2 ' 
j dAE A (A)F A (<D 



(1) 



where Aa is the extinction in the A band, E\ is the effective area 
of that band, and F,\ is the source flux density. The result is a 
Galactic extinction of 1.10, 1.44, 1.74, 2.40, 3.04, and 2.92 mag 
from the v to the uvwl band, respectively. 

We verified the stability of our results by iterating the pro- 
cedure with the recalibrated UVOT flux densities and A e $. The 
bottom panel of Fig.|4]shows the mean SED obtained after recal- 
ibration of the UVOT data according to our procedure. For com- 
parison, we also show the OM and UVOT SEDs derived from the 
average spectra shown in the top panel, for which the amount of 
Galactic extinction was calculated from the lCardelli et al.l |l989) 
law at the standard A e ff. Notice that the recalibration process has 
shifted Aes (uvwl) redward so much that it overlaps with the stan- 
dard A e s(uv ml). This is becau se of the noticeable red tail of the 
uvwl filter dPoole et all 120081) as well as to the red BL Lacertae 
spectrum. Recalibration has solved the uvwl -dip problem, which 
is a common feature of UVO T SEDs for a number of blazars 
at different redshifts (see e.g. |Villata et all l2008t iRaiteri et all 
HooUd! Ammando et a n i2oo9ir Moreover, it seems to confirm 

the UV excess claimed bv lRaiteri et al.l (120091) that was ascribed 
to thermal emission from the accretion disc, even if this excess 
may be less pronounced than indicated by the OM data. Our 
analysis highlights the importance of calculating the amount of 
extinction in the critical UV bands, close to the 2175 A bump, 
by folding the Galactic mean extinction law through the effec- 
tiv e area curves and sou r ce sp ectrum. In any case, as pointed out 
bv lFitzpatrick & Massa (2007), one has to keep in mind that the 
use of an average dereddening curve implies a significant error 
owing to the scatter of Galactic extinction curves. 



6 N otice that our CF(wvwl) corresponds to that given by Poole et al. 
(2008) for the Pickles stars. 
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Fig. 5. Swift-XRT spectrum of BL Lacertae on 2008 August 29. 
The bottom panel shows the deviations of the observed data from 
the folded model (a single power law with fixed absorption) in 
units of standard deviations. 



3.2. XRTdata 



The X-ray Telescope (XRT; iBurrows et al . 2005) data were pro- 
cessed with version 0.12.3 of the xrtpipeline task contained 
in the FTOOLS package, applying standard screening criteria. 
Inspection of the light curves revealed that the count rate was 
low, from 0. 15 to 0.24 counts s _1 , so that observations were per- 
formed in photon counting mode, and no correction for pile- 
up was necessary. Source and background spectra were ex- 
tracted with xselect from a circular region of 20 pixel (47 
arcsec) radius centred on the source and from a surrounding 
annulus of 30 and 50 pixel radii, respectively. We used ver- 
sion 011 of the response matrix available in the HEASARC 
calibration database (CALDB), and calculated the ancillary re- 
sponse file with xrtmkarf, using the exposure map created by 
xrtpipeline. The source spectra were binned with grppha to 
have a minimum of 20 counts in each bin, and they were finally 
analysed with version 12.5.1 of the Xspec task, using the energy 
channels greater than 0.3 keV. 

Spe ctral analysis was p erformed for each observation fol- 
lowing l^iterr^raL J200^: we first fitted a single power law 
with free absorptiory, and then fixed the Galactic absorption to 
Nh = 3.4 x 10 21 ctrT 2 , which takes into account both atomic and 
molecular column density. Statistics is not good enough to eval- 
uate if a double power law model can improve the fit. The results 
of spectral fitting on XRT data are reported in Table|2]for all ob- 
servations with an exposure longer than 3 ksec; Col. 1 gives the 
date and start time of the observation; Col. 2 its duration; Col. 3 
the hydrogen column; Col. 4 the power law photon index; Col. 5 
the 1 keV flux density; Col. 6 the^- 2 /v (and degrees of freedom). 
One spectrum (August 29) is shown in Fig. [5] 

Fits with free absorption resulted in a very variable A^h, 
which is unlikely to correspond to a real change of absorption. 
The average and median values are 3.46 and 3.44 x 10 21 cm -2 , 
respectively, confirming that the value assumed for the Galactic 
absorption is quite reasonable. Hence, we favoured the second 
model, whose y 1 l v is usually smaller than in the A^-free case, 
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Fig. 6. Details of the multiwavelength behaviour of BL Lacertae 
in the period around the Swift observations. The X-ray flux den- 
sity (at 1 keV, top) is compared to that in the UV (uvwl band), 
optical (R band), near-IR (K band), and radio (37 GHz, bottom) 
frequency range. 



We adopted the Tuebingen-Boulder ISM absorption model 



1 We ado pted tt 
dWilms et aUl200Ch . 



and that produces results with smaller errors (because of one 
degree of freedom more). In only two cases (August 23 and 
September 2) a double power law model with absorption fixed 
to the Galactic value clearly improved the fit. 

The photon index T ranges from 1.92 to 2.25, indicating 
a spectrum that oscillates from moderately hard to moderately 
soft. The average value is 2.07, with standard deviation of 0.08. 
To understand whether these spectral changes correspond to real 
variations or are owing to noise, we recall the definition of the 
mean fractional variation F val = Vcr 2 - <5 2 / < / > dPetersonl 
2001), which is commonly used to characterise variability. Here 
< / > is the mean value of the variable we are analysing, cr 2 
its variance, and 6 1 the mean square uncertainty. In our case, 
cr 2 = 0.006 is smaller than 6 2 = 0.01 1, so that the result is imag- 
inary; thus we conclude that the variations are consistent with 
noise rather than source variability. 

The 1 keV flux density varies between 1.34 and 2.24 yuJy, 
with a mean value of 1 .75 and standard deviation of 0.24. In this 
case F v - dr =0.11, and the variations can be considered reliable. 

Multiwavelength light curves of BL Lacertae in the period 
around the Swift observations are shown in Fig. [6] The source 
behaviour at 1 ke V differs from the common trend characterising 
the UV, optical, and near-IR bands. In particular, the X-ray flux 
peaks when the near-IR-UV fluxes reach a minimum. However, 
there are also similarities, like the flux increase at the beginning 
of the common observing period, and the final decrease. This 
may indicate that the 1 keV flux behaviour sometimes is related 
to the brightness changes that occur at lower wavelengths, while 
in other cas es another va riability mechanism prevails. Indeed, 
according to lRaiteri et al.l(l2009l) this frequency domain receives 
the variable contribution of two different emission components 
(see also Sect. [5J. 
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Table 2. Spectral fitting to the XRT data from Swift observations in 2008. Only exposures longer than 3 ksec are considered. For 
each epoch the first line reports the result of the single power law model with free absorption, while the second line shows that 
obtained when fixing the Galactic absorption to Nu = 3.4 x 10 21 cirT 2 . 
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Fig. 7. UV and X-ray light curves obtained by the OM and EPIC 
instruments onboard XMM-Newton on 2008 May 16-17. Black, 
red, and blue filled circles represent UW1, UVM2, and UVW2 
magnitudes, respectively. Arrows mark the data points used for 
the SED shown in Fig. [9] Green diamonds, black squares, and 
red triangles refer to pn (0.5-10 keV), MOS1 (0.3-10 keV), and 
MOS2 (0.3-10 keV) count rates, respectively. Plus signs show 
the corresponding backgrounds. In all panels horizontal dashed 
lines indicate average values and dotted lines represent standard 
deviations from the mean. 



4. XMM-Newton observations 

The X-ray Multi-Mirror Mission (XMM) - Newton satellite ob- 
served the source during revolution 1545, on 2008 May 16-17, 
with a total exposure of ~ 134ks. Data were processed with the 
Science Analysis System (SAS) package version 9.0. 



4.1. OMdata 



The Optical Monitor (OM: lMason et al.ll200ll) onboard XMM- 
Newton is a 30-cm telescope carrying six optical/UV filters, and 
two grisms. BL Lacertae observations in May 2008 consisted of 
10 subsequent exposures in UVWl, followed by 9 in UVM2, 
and then 8 in UVW2. All exposures were ~ 4000 s long. We 
used the SAS task omi chain to reduce the data and the tasks 
omsource and omphotom to derive the source magnitude. The 
error on the aperture photometry is 0.03, 0.04, and 0.09 mag 
for the UVWl, UVM2, and UVW2 filters, respectively. The re- 
sulting light curves are shown in Fig.|7j average magnitudes are 
UVWl = 15.45, UVM2 = 16.25, and UVW2 = 16.54. 

To obtain flux densities for further analysis, OM magnitudes 
were c orrected for the Gala ctic extinction calculated according 
to the ICardelli et all (1 19891 see Sect. 13 laws at the effective 
wavelengths of the OM filters (2910, 2310, and 2120 A for the 
UVWl, UVM2, and UVW2 filters, respectively). Conversion of 





Win r 



0.5 1 2 5 10 

Energy (keV) 

Fig. 8. EPIC spectrum of BL Lacertae on 2008 May 16-17; 
black squares, red triangles, and green diamonds represent 
MOS1, MOS2, and pn data, respectively. The bottom panel 
shows the deviations of the observed data from the folded model 
(a double power law with fixed absorption) in units of standard 
deviations. 



de-reddened magnitudes into flux densities was done with re- 
spect to Vega. 

4.2. EPIC data 

The European Photon Imaging Camera (EPIC) on board XMM- 
Newt on includ es three detectors: M OS1, MOS2 dTurner et all 
1200 lb . and pn (Stri ider et al.L l200"lb . All instruments were used 
with a thin filter. The two MOS cameras observed in small- 
window imaging mode, while pn was used in timing mode. 

We followed the standard prescription to reduce the data, in- 
cluding filtering of high background periods with a threshold of 
0.35 counts s -1 for MOS, but with a stricter threshold of 0.1 
counts s" 1 forpn. 

For both MOS1 and MOS2, we created a filtered sky image, 
and extracted the source counts from a 50 arcsec radius circular 
region, while background was evaluated in a circle on an ex- 
ternal CCD. As for pn, we extracted the source counts from a 
strip between RAWX=35 and 39, and the background from two 
strips at columns 24-28 and 48-52. To get the most reliable and 
best calibrated events, we used the FLAG==0 selection expres- 
sion and kept only single and double events (PATTERN < =4). 
We verified that pile-up effects were not affecting the MOS data 
with the epatplot task. 

Through the grppha task of the FTOOL package we binned 
the source spectra to have a minimum of 25 counts in each bin 
and then analysed them together by means of the Xspec task 
of the XANADU package. Only spectral bins corresponding to 
energies between 0.3 and 10 keV for MOS1 and MOS2 and in 
the range 0.5-10 keV for pn were considered, because they have 
both a better calibration and a higher signal-to-noise ratio. 

The three EPIC spectra were analysed together by first fit- 
ting a single power law with free absorption, and then fixing the 
Galactic absorption to Nn = 3.4 x 10 21 crrT 2 . We also tried a 
double power law with the same Galactic absorption. The re- 
sults are shown in Table [3] (see also Fig. [8}. In agreement with 
Raite ri et al . (2009), the^ 2 /v suggests that a single power law 
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Table 3. Spectral fitting to the EPIC (pn+MOS 1 +MOS2) data 
from the XMM-Newton observations of 2008 May 16-17. The 
first line reports the result of the single power law model with 
free absorption; the second shows that obtained when fixing the 
Galactic absorption to Nn = 3.4 X 10 21 cm" 2 , and the third that 
from a double power law model with the same Galactic absorp- 
tion. 



Table 4. Main parameters of the helical jet model for the fit to 
both the 2008 faint-state SED and the 1997 outburst- state SED. 
For each epoch, "low" and "high" refer to the low- and high- 
energy synchrotron plus self inverse-Compton components, re- 
spectively. 
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with Galactic absorption of Nn = 3.4 x 10 21 crrT 2 does not rep- 
resent a good fit to the data. As for the other two fits, the double 
power law seems to better fit the data, which is also confirmed 
by a very low F-test probability of ~ 6.8 x 10~ 16 . This implies a 
strong spectral curvature. 

To check for possible flux variations, we extracted X-ray 
light curves from the same source and background regions de- 
fined for the spectra, with the same selection expressions. We 
considered only the events in the time intervals free from high 
background and belonging to the 0.3-10 keV energy range for 
MOS1 and MOS2, and 0.5-10 keV for pn. The source counts 
were corrected for the background and then binned in one hour 
intervals. The results are shown in Fig. [7] which also displays 
the behaviour of the background to check the reliability of the 
flux variations. The background increased significantly only in 
the last 6 hours. Just before, at JD ~ 2454604.05, there is a 
small flare clearly visible in all three light curves. Mean source 
rates for the whole period are 1.63, 0.61, and 0.64 counts s~' for 
pn, MOS1, and MOS2, respectively, with standard deviations of 
0.15, 0.04, and 0.03 counts s" 1 . This means fractional variations 
,F V ar of 7%, 6%, and 4%, which do not change significantly if we 
exclude the last 6 hours. Hence, we can conclude that the X-ray 
flux of BL Lacertae is mildly variable on an hour time scale. 



5. Modelling the SED 

Figure [9] shows the broad-band SED of BL Lacertae in differ- 
ent brightness states. The SED corresponding to 2008 May 16- 
17 includes the XMM-Newton UV and X-ray data analysed in 
Sect. |4] In order to avoid offsets caused by source variability, the 
OM spectrum was constructed with the last UVW1 data point, 
the first UVM2 (which is close to the last UVM2), and the first 
UVW2 datum of the corresponding light curves (see Fig. [7]). 
These data indicate a hard UV spectrum. Two SEDs in the figure 
refer to the Swift observations of 2008 August 25 and 30, which 
were chosen among those analysed in Sect.[3]because of the dif- 
ferent spectral slope in the X-ray band. The low-frequency part 
of these three SEDs is built with GASP optical and radio data 
(see Sect. |2J; for August 30 near-IR data were also available. 
The strong source variability in the optical band requires that the 
optical (and near-IR) data are simultaneous with the satellite ob- 
servations. In the radio bands flux variations are slower, so that 
we used data taken within 2-3 days from the satellite ones, when 
simultaneous data were not available. Near-IR, optical, and UV 
data were corrected for the Galactic extinction; the near-IR and 
optical data were also corrected for the contribution of the host 
galaxy (see Sect. [2]). In August 2008 the Fermi y-ray satellite 
detected BL Lacertae; th e Fermi data we plotted in Fig. [9] were 
derived from Abd o et al 1 (l2010ah . 
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The August 2008 SEDs indicate a faint, synchrotron- 
dominated state o f the source that we fitted with th e rotating he- 
lical jet model by Villata & RaiterU ( T9991 see also lRaiteri et al] 
1 19991 iRaiteri et al.ll2003l lOstorero et al.ll2004l) . This model has 
been used by IRaiteri et al.l d2009l) to fit the broad-band SED of 
BL Lacertae in December 2007 - January 2008. Their main find- 
ing was that the BL Lacertae broad-band SED cannot be ex- 
plained by a single synchrotron component plus its self inverse- 
Compton emission. Indeed, the very strong historical X-ray vari- 
ability requires an additional synchrotron (plus self inverse- 
Compton) component. Moreover, the UV excess suggests ther- 
mal contribution from the accretion disc (see Sect. 13. U . The 
SED analysed by Raiteri TtaD d2009t) lacked simultaneous y-ray 
data, which made it impossible to constrain the inverse-Compton 
emission of the high-energy component. The 2008 August SED 
in Fig.|9]now offers us the possibility to perform a more detailed 
analysis. 

In addition, we also display in Fig. [9] the broad-band SED 
corresponding to the big outburst of July 1997, which showed 
a considerable inverse-Compton dominance. The X-ray spectra 
plotted in the figure are the res ult of the com b ined a nalysis of 
the ASCA and RXTE data by iTanihata et ail j2000h . Because 
of the very strong variability of the source in that period, the 
authors distinguished between a low state, which was well fit- 
ted by a single power law model, and a flare state, for which 
the best fit was obtained with a double power law model. This 
last fit resulted in a very strong spectral curvatur^. In July 1997 
observations in the y-ray band were performed by CGRO. The 
data from the EGR ET instrum ent onboard CGRO in Fig. [9] were 
taken from Bloom et al. (1997), while those from the OSSE de- 
tector were derived from the High Energy Astrophysics Science 
Archive Research Centei0 (HEASARC). The low-frequency in- 
formation is from the WEBT archive; the range of optical flux 
variation in the period is indicated. This outburst state of the 
source was fitted with the same rotating helical model that we 
used to fit the faint state of 2008. 



8 We notice that Tanihataet al. (2000) adopted a Galactic total ab- 
sorption of N H = 4.6 x 10 21 cm~ 2 . 

9 http://heasarc.gsfc.nasa.gov/ 
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Fig. 9. Broad-band SEDs of BL Lacertae in 
August 2008 (blue) and July 1997 (red). The 
2008 SED is built with UV and X-ray data 
from two epochs of Swift observations (Sect. 
0, 2008 August 25 and 30, together with 
simultaneous optical-to-radio data from the 
GASP-WEBT collaboration, and contempora- 
neous y-ray data from Fermi (from lAbdo et all 
I2010ah. The 1997 S EP in cludes EGRET data 
from Bloom et al.1 dl997t> . OSSE data from 
the HEA SARC archive ASCA +RXTE spec- 
tra from iTanihata et a l. (2000), while low- 
frequency data are from the WEBT archive. 
Solid lines represent model fits; we distin- 
guish the low-energy emission component (dot- 
ted lines) from the high-energy one (dashed 
lines); the contribution by an accretion disc of 
~ 16000 K and 5 x 10 44 erg s -1 is marked with 
a dotted-dashed line. We also show in black the 
SED corresponding to the XMM-Newton ob- 
servations of May 2008 (Sect. [4]l; both the sin- 
gle power law with free absorption and double 
power law with Galactic absorption fits to the 
EPIC spectra are displayed. 



In performing the model fits our aim was to see whether it 
was possible to reproduce the high and low states by changing 
the geometrical con figuration only. Moreo ver, we took into ac- 
count the results by lLarionov et al.l d2010i) . who explained the 
long-term BL Lacertae optical and near-IR variability in terms 
of variations of the Doppler boosting factor due to changes of 
the viewing angle of the emitting region. 

The resulting model parameters are reported in Table |4] 
while the corresponding fits are shown in Fig. [9] We also in- 
cluded blackbody radiation from an accretion disc with a lu- 
minosity of 5 x 10 44 ergs _1 and a temperature of ~ 16000 K. 
The low-energy (radio-to-optical and related inverse-Compton) 
emission component comes from a helix portion with a pitch an- 
gle £ = 2°, covering an angle a = 180°. The maximum Lorentz 
factor of the relativistic electrons is log y max (0) = 3.7, while the 
bulk Lorentz factor of the plasma in the jet is T = 7. The high- 
energy (UV-X-ray and related inverse-Compton) emission com- 
ponent comes from a helix portion with a pitch angle £ = 8°, 
covering an angle a = 360°. The maximum Lorentz factor of the 
relativistic electrons is logy max (0) = 3.8, while the bulk Lorentz 
factor is F = 13. These parameters, as well as those defining 
the power laws according to which the maximum and minimum 
emitted frequencies and the flux densities decrease with distance 
from the jet apex, are maintained fixed. The difference between 
the fits to the 1997 and 2008 SEDs is only due to a variation of 
the orientation of the jet emitting regions, through only two geo- 
metric parameters: the angle between the helix axis and the line 
of sight if/ and the rotation angle 0. The outburst state requires 
a better alignment of the emitting regions with the line of sight, 
which implies if/ approaching the helix pitch angle, and smaller 
rotation angles (f>. 

6. Discussion 

The satisfactory fits that we obtained in the previous section for 
both the outburst and faint states of BL Lacertae give strength 
to our interpretation of the source SED in terms of two syn- 
chrotron plus self inverse-Compton emission components. In the 
SED, the synchrotron peak of the low-energy component falls in 
the near-IR band, and its inverse-Compton reaches a maximum 




Fig. 10. Sketch of our helical jet model during both the 2008 
faint state and the 1997 outburst state. The angle between the jet 
axis and the line of sight, if/, has been multiplied by a factor 10 
with respect to the values given in Table|4]for clarity. We distin- 
guish the inner region, emitting the high-energy synchrotron plus 
self inverse-Compton component (purple-blue) from the outer 
zone, where the low-energy emission component is produced 
(green-yellow -red). 



in the 1-50 MeV energy range. The synchrotron and inverse- 
Compton peaks of the high-energy component occur in the far- 
UV-soft-X-ray band and in the energy range 0.2-5 GeV, respec- 
tively. Whether these two components come from two distinct 
helices or from different regions inside the same helical jet is not 
clear. We consider it more likely that there is a unique jet, where 
the high-energy component comes from a region closer to the 
emitting jet apex than the low-energy one. We notice that a dou- 
ble synchrotron component is not an unusual interpretation for 
the blazar emission properties , as it has been propos ed also for 
Mkn 421 dDonnarumma et all [2009b and 3C 454.3 dOrie et all 
120101) . 
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Moreover, the fits show that the whole range of BL Lacertae 
multiwavelength variability can be interpreted in terms of orien- 
tation effects. Although the rotating helical jet model we have 
adopted in the previous section is not a physically complete 
model, but more a phenomenological approach, it has the ad- 
vantage of taking into account variations of the orientation of 
the emitting regions with respect to the line of sight, with conse- 
quent changes of the Doppler beaming factor. This is an aspect 
that is usually neglected by theoretical models of blazar emis- 
sion, which explain flux and spectral changes uniquely in terms 
of energetic processes inside the jet. 

Our interpretatio n is in line with previous results. 
Mars cher et al.1 (120081) analysed the evolution of the BL Lacertae 
optical polarization during 2005, and suggested th at the plasma 
flows a long helical streamlines. According to IVillata et al.l 
(2009b), the long-term optical and radio behaviour of BL 
Lacertae suggests a scenario where the emitting plasma flows 
along a rotating helical path in a curved jet. This rotating he- 
lical structure could be caused by orbital motion in a binary 
black hole system, coupled with the interaction of the plasma 
jet with the surrounding medium. Indeed, the binary black hole 
scenario could explain th e periodicity observed in th e radio light 
curves of BL Lacertae (IVillata et al.Ll2004al l2009b). the discov- 
eryof a precessing jet nozzle with the VLBA (Stirling et al., 
2003), a nd possibly t he pa rsec-to-kiloparsec jet misalignment 
(see e.g. iKharb et all 120101) . Moreover, the ana lysis of the BL 
Lacert ae spectral evolution in 2000-2008 by lLarionov et alj 
(2010) favoured a picture where the optical and near-IR flux and 
colour variability can be explained by a variable viewing angle 
of the emitting region. These authors also suggested that a fractal 
helical structure may be at the origin of the different time scales 
of variability. 

The values of the angles £ = 2° and if/ - 4.5-5°, as well as 
the Lorentz factor F — 7 adopted in the model fits for the low- 
energ y component, agree ve ry well with the corresponding val- 
ues of lLarionov et al. d2010l) . thus supporting the common inter- 
pretation. The helix pitch angle £ = 8° found for the high-energy 
component indicates that this inner jet helical region would be 
more twisted than the outer, lower-energy one. In practice, our 
model results indicate a helical jet whose axis is bent between 
the X-ray and optical regions by about 2° (see the values of if/ in 
Table |4]i and that is more wrapped near the apex and then tends 
to relax with a decreasing pitch angle. The different orientations 
assumed by such a jet in 2008 and 1997 are sketched in Fig.fTUl 
where all the angles are strongly increased for clarity. 

The thermal emission component that we added to the he- 
lical jet model is justified by the UV excess found in the OM 
data from XMM-Newton, and (though with less evidence) in the 
UVOT data from Sw ift. As discussed in the present paper and in 
Raiteri et al. (2009), the amount of this excess strongly depends 
on both the Galactic extinction and instrument cali bration, but it 
is not easy to cancel it out completely. In any case, Capetti et al. 
(2010) showed that after twelve y ears from the fir s t det ection 
of the Ha broad em ission li ne by Vermeulen et al.l d 1995b and 
ICorbett et all (1 1 996l see also lCorbett et al.ll2000h ."the Ha line is 
still there, even more luminous than before. This suggests that a 
disc is also there to photoionise the broad line region. 

Photons coming from the disc or broad line region could then 
enter the jet, and be inverse-Compton scattered, giving rise to 
other high-energy emission components that are sometimes in- 
voked to account for the SED properties of blazars. In particu- 
lar, the 1997 o utburst state has previously been interpreted by 
Madejski etaD (1 1999b in terms of three emission components: 
synchrotron, synchrotron self-Compton, and Comptonisation of 



the broad emission line fl ux. Simil ar results were obtain ed by 
iBottcher & Blooml (l2000l) and by iRavasio et all (120021) . Our 
"geometrical" interpretation does not require these external- 
Compton emission components, which are not expected to con- 
tribute if the jet emission regions are parsecs away from the 
central black hole (see e .g. ISikora et all 120081; FMarscher et al.L 
l2010HAbdo eTafll2010bh . 
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